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(54) Method and system for recognizing, indexing, and searching acoustic signals 



(57) A computerized method extracts features from 
an acoustic signal generated from one or more sources. 
The acoustic signal are first windowed and filtered to 
produce a spectral envelope for each source. The di- 
mensionality of the spectral envelope is then reduced to 
produce a set of features for the acoustic signal. The 
features in the set are clustered to produce a group of 
features for each of the sources. The features in each 



group include spectral features and corresponding tem- 
poral features characterizing each source. Each group 
of features is a quantitative descriptor that is also asso- 
ciated with a qualitative descriptor. Hidden Markov mod- 
els are trained with sets of known features and stored 
in a database. The database can then be indexed by 
sets of unknown features to select or recognize like 
acoustic signals. 
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Description 

RELATED APPLICATION 

5 [0001] This application is a Continuation-in-Part Application of U.S. Patent Application 09/346,854 "Method for Ex- 
tracting Features from a Mixture of signals" filed by Casey on July 2, 1999. 

FIELD OF THE INVENTION 

10 [0002] The invention relates generally to the field of acoustic signal processing, and in particular to recognizing, 
indexing and searching acoustic signals. 

BACKGROUND OF THE INVENTION 

15 [0003] To date, very little work has been done on characterizing environmental and ambient sounds. Most prior art 
acoustic signal representation methods have focused on human speech and music. However, there are no good rep- 
resentation methods for many sound effects heard in films, television, video games, and virtual environments, such 
footsteps, traffic, doors slamming, laser guns, hammering, smashing, thunder claps, leaves rustling, water spilling, etc. 
These environmental acoustic signals are generally much harder to characterize than speech and music because they 

20 often comprise multiple noisy and textured components, as well as higher-order structural components such as itera- 
tions and scattering. 

[0004] One particular application that could use such a representation scheme is video processing. Methods are 

available for extracting, compressing, searching, and classifying video objects, see for example the various MPEG 
standards. No such methods exist for "audio" objects, other than when the audio objects are speech. For example, it 
25 maybe desired to search through a video library to locate all video segments where John Wayne is galloping on a 
horse while firing his six-shooter. Certainly it is possible to visually identify John Wayne or a horse. But it much more 
difficult to pick out the rhythmic clippidy-clop of a galloping horse, and the staccato percussion of a revolver. Recognition 
of audio events can delineate action in video. 

[0005] Another application that could use the representation is sound synthesis. It is not until the features of a sound 
30 are identified before it becomes possible to synthetically generate a sound, other than be trial and error. 

[0006] In the prior art, representations for non-speech sounds have usually focused on particular classes of non- 
speech sound, for example, simulating and identifying specific musical instruments, distinguishing submarine sounds 
from ambient sea sounds and recognition of underwater mammals by their utterances. Each of these applications 
requires a particular arrangement of acoustic features that do not generalize beyond the specific application. 
35 [0007] In addition to these specific applications, other work has focused on developing generalized acoustic scene 
analysis representations. This research has become known as "computational auditory scene analysis." These systems 
require a lot of computational effort due to their algorithmic complexity. Typically, they use heuristic schemes from 
Artificial Intelligence as well as various inference schemes. 

[0008] Whilst such systems provide valuable insight into the difficult problem of acoustic representations, the per- 
40 formance of such systems has never been demonstrated to be satisfactory with respect to classification and synthesis 
of acoustic signals in a mixture. 

[0009] In yet another application, sound representations could be used to index audio media including a wide range 
of sound phenomena including environmental sounds, background noises, sound effects (Foley sounds), animal 
sounds, speech, non-speech utterances and music. This would allow one to design sound recognition tools for search- 
es ing audio media using automatically extracted indexes. Using these tools, rich sound tracks, such as films or news 
programs, could be searched by semantic descriptions of content or by similarity to a target audio query. For example, 
it is desired to locate all film clips where lions roar, or elephants trumpet. 

[0010] There are many possible approaches to automatic classification and indexing. Wold et al.," IEEE Multimedia, 
pp. 27-36, 1996, Martin et al., "Musical instrument identification: a pattern-recognition approacin," Presented at the 

50 136th Meeting of the Acoustical Society of America, Norfolk, VA, 1998, describe classification strictly for musical in- 
struments. Zhang et al., "Content-based classification and retrieval of audio," SP\B 43''^ Annual Meeting, Conference 
on Advanced Signal Processing Algorithms, Architectures and Implementations VIII, 1998, describes a system that 
trains models with spectrogram data, and Boreczky et al., "A hidden Markov model framework for video segmentation 
using audio and image feafures," Proceedings of ICASSP'98, pp.3741 -3744, 1998 uses Markov models. 

55 [0011] Indexing and searching audio media is particularly germane to the newly emerging MPEG-7 standard for 
multimedia. The standard needs a unified interface for general sound classes. Encoder compatibility is a factor in the 
design. Then, a "sound" database with indexes provided by one implementation could be compared with those extracted 
by a different implementation. 
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SUMMARY OF THE INVENTION 

[0012] A computerized method extracts features from an acoustic signal generated from one or more sources. The 
acoustic signal are first windowed and filtered to produce a spectral envelope for each source. The dimensionality of 

the spectral envelope is then reduced to produce a set of features for the acoustic signal. The features in the set are 
clustered to produce a group of features for each of the sources. The features in each group include spectral features 
and corresponding temporal features characterizing each source. 

[001 3] Each group of features is a quantitative descriptor that is also associated with a qualitative descriptor. Hidden 
Markov models are trained with sets of known features and stored in a database. The database can then be indexed 
by sets of unknown features to select or recognize like acoustic signals. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] 

Figure 1 is a flow diagram of a method for extracting features from a mixture of signals according to the invention; 

Figure 2 is a block diagram of the filtering and windowing steps; 

Figure 3 is a block diagram of normalizing, reducing , and extracting steps; 

Figures 4 and 5 are graphs of features of a metallic shaker; 

Figure 6 is a block diagram of a description model for dogs barking; 

Figure 7 is a block diagram of a description model for pet sounds; 

Figure 8 is a spectrogram reconstructed from four spectral basis functions and basis projections; 

Figure 9a is a basis projection envelope for laughter; 

Figure 9b is an audio spectrum for the laughter of Figure 9 

Figure 10a is a log scale spectrogram for laughter; 

Figure 10b is a reconstructed spectrogram for laughter; 

Figure 11a is a log spectrogram for dog barking; 

Figure 1 1 b is a sound model state path sequence of states through a continuous hidden Markov model for the dog 
barking of Figure 11a; 

Figure 12 is a block diagram of a sound recognition classifier; 

Figure 13 is a block diagram of a system for extracting sounds according to the invention; 

Figure 14 is a block diagram of a process for training a hidden Markov model according to the invention; 

Figure 15 is a block diagram of a system for identifying and classifying sounds according to the invention; 

Figure 1 6 is a graph of a performance of the system of Figure 1 5; 

Figure 17 is a block diagram of a sound query system according to the invention; 

Figure 1 8a is a block diagram of a state path of laugher; 



Figure 18b is a state path histograms of laughter; 
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Figure 19a are state paths of matching laughters; and 

Figure 19b are state path histograms of matching laughters. 

5 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0015] Figure 1 shows a method 1 00 for extracting spectral and temporal features 1 08-1 09 from a mixture of signals 
101 according to my invention. My method 100 can be used for characterizing and extracting features from sound 
recordings for classification of the sound sources and for re-purposing in structured multi-media applications such as 
10 parametric synthesis. The method can also be used to extract features from other linear mixtures, or for that matter 
from multi-dimensional mixtures. The mixture can be obtained from a single source, or from multiple sources such as 
a stereo sound source. 

[0016] In order to extract features from recorded signals, I use statistical techniques based on independent compo- 
nent analysis (ICA). Using a contrast function defined on cumulative expansions up to a fourth order, the ICA transform 
15 generates a rotation of the basis of the time-frequency observation matrices 1 21 . 

[0017] The resulting basis components are as statistically independent as possible and characterize the structure 
of the individual features, e.g., sounds, within the mixture source 101 . These characteristic structures can be used to 
classify the signal, or to specify new signals with predictable features. 

[0018] The representation according to my invention is capable of synthesizing multiple sound behaviors from a 
20 small set of features. It is able to synthesize complex acoustic event structures such as impacts, bounces, smashes 
and scraping as well as acoustic object properties such as materials, size and shape. 

[0019] In the method 100, the audio mixture 101 is first processed by a bank of logarithmic filters 110. Each of the 
filters produces a band-pass signal 111 for a predetermined frequency range. Typically, forty to fifty band-pass signals 
111 are produced with more signals at lower frequency ranges than higher frequency ranges to mimic the frequency 
25 response characteristics of the human ear. Alternatively, the filters can be a constant-Q (CQ) or wavelet filterbank, or 

they can be linearly spaced as in a short time fast Fourier transform representation (STFT). 

[0020] In step 120, each of the band-pass signals is "windowed" into short, for example, 20 millisecond segments 
to produce observation matrices. Each matrix can include hundreds of samples. The details of steps 110 and 120 are 
shown in greater detail in Figures 2 and 3. It should be noted that the windowing can be done before the filtering. 
30 [0021] In step 130 a singular value decomposition (SVD) is applied to the observation matrices 121 to produce 
reduced dimensionality of the matrices 131. SVD were first described by the Italian geometer Beltrami in 1873. The 
singular value decomposition is a well-defined generalization of the principal component analysis (RCA). The singular 
value decomposition of an m x r? matrix is any factorization of the form: 

35 

x=ui: Y 

where is an m x m orthogonal matrix; i.e. L/has orthonormal columns, \^is an n xn orthogonal matrix, and E is an 
m X n diagonal matrix of singular values with components Oy = 0 if / is not equal to / 

40 [0022] As an advantage and in contrast with PGA, the SVD can decomposes a non-square matrix, thus it is possible 
to directly decompose the observation matrices in either spectral or temporal orientation without the need for a calcu- 
lating a covariance matrix. Because the SVD decomposes a non-square matrix directly, without the need for a covar- 
iance matrix, the resulting basis is not as susceptible to dynamic range problems as the PGA. 
[0023] I apply an optional independent component analysis (IGA) in step 1 40 to the reduced dimensionality matrices 

45 131 . An IGA that uses an iterative on-line algorithm based on a neuro-mimetic architecture for blind signal separation 
is well known. Recently, many neural- network architectures have been proposed for solving the IGA problem, see for 
example, U.S. Patent No. 5,383,1 64 "Adaptive system for broadband multisignal discrimination in a channel with re- 
verberation," issued to Sejnowski on January 17, 1995. 

[0024] The IGA produces the spectral and temporal features 1 08-1 09. The spectral features, expressed as vectors, 
50 correspond to estimates of the statistically most independent component within a segmentation window. The temporal 
features, also expressed as vectors, described the evolution of the spectral components during the course of the 
segment. 

[0025] Each pair of spectral and temporal vectors can be combined using a vector outer product to reconstruct a 
partial spectrum for the given input spectrum. If these spectra are invertable, as a filterbank representation would be, 
55 then the independent time-domain signals can be estimated. For each of the independent components described in 
the scheme, a matrix of compatibility scores for components in the prior segment is made available. This allows tracking 
of components through time by estimating the most likely successive correspondences. Identical to the backward 
compatibility matrix, only looking forward in time. 
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[0026] An independent components decomposition of an audio track can be used to estimate individual signal com- 
ponents within an audio track. Whilst the separation problem is intractable unless a full-rank signal matrix is available 
(N linear mixes of N sources), the use of independent components of short temporal sections of frequency-domain 
representations can give approximations to the underlying sources. These approximations can be used for classification 
5 and recognition tasks, as well as comparisons between sounds. 

[0027] As shown in Figure 3, the time frequency distribution (TFD) can be normalized by the power spectral density 
(PSD) 115 to diminish the contribution of lower frequency components that carry more energy in some acoustic do- 
mains. 

[0028] Figures 4 and 5 respectively show the temporal and spatial decomposition for a percussion shaker instrument 

10 played at a regular rhythm. The observable structures reveal wide-band articulate components corresponding to the 
shakes, and horizontal stratification corresponding to the ringing of the metal shell. 

Applications for Acoustic Features of Sounds 

15 [0029] My invention can be used in a number of applications. The extracted features can be considered as separable 
components of an acoustic mixture representing the inherent structure within the source mixture. Extracted features 
can be compared against a set of a-priori classes, determined by pattern-recognition techniques, in order to recognize 
or identify the components. These classifiers can be in the domain of speech phonemes, sound effects, musical in- 
struments, animal sounds or any other corpus-based analytic models. Extracted features can be re-synthesized inde- 

20 pendently using an inverse filter-bank thus achieving an "unmixing" of the source acoustic mixture. An example use 
separates the singer, drums and guitars from an acoustic recording in order to re-purpose some components or to 
automatically analyze the musical structure. Another example separates an actor's voice from background noise in 
order to pass the cleaned speech signal to a speech recognizer for automatic transcription of a movie. 
[0030] The spectral features and temporal features can be considered separately in order to identify various prop- 

25 erties of the acoustic structure of individual sound objects within a mixture. Spectral features can delineate such prop- 
erties are materials, size, shape whereas temporal features can delineate behaviors such as bouncing, breaking and 
smashing. Thus a glass smashing can be distinguished from a glass bouncing, or a clay pot smashing. Extracted 
features can be altered and re-synthesized in order to produce modified synthetic instances of the source sound. If 
the input sound is a single sound event comprising a plurality of acoustic features, such as a glass smash, then the 

30 individual features can be controlled for resynthesis. This is useful for model-based media applications such as gen- 
erating sound in virtual environments. 

Indexing and Searching 

35 [0031] My invention can also be used to index and search a large multimedia database including many different 
types of sounds, e.g., sound effects, animal sounds, musical instruments, voices, textures, environmental sounds, 
male sounds, female sounds. 

[0032] In this context, sound descriptions are generally divided into two types: qualitative text-based description by 

category labels, and quantitative description using probabilistic model states. Category labels provide qualitative in- 
40 formation about sound content. Descriptions in this form are suitable for text-based query applications, such as Internet 
search engines, or any processing tool that uses text fields. 

[0033] In contrast, the quantitative descriptors include a compact information about an audio segment and can be 
used for numerical evaluation of sound similarity. For example, these descriptors can be used to identify specific in- 
struments in a video or audio recording. The qualitative and quantitative descriptors are well suited to audio query-by- 
45 example search applications. 

Sound Recognition Descriptors and Description Schemes 

Qualitative Descriptors 

50 

[0034] While segmenting an audio recording into classes, it is desired to gain pertinent semantic information about 

the content. For example, recognizing a scream in a video soundtrack can indicate horror or danger, and laughter can 
indicate comedy. Furthermore, sounds can indicate the presence of a person and therefore the video segments to 
which these sounds belong can be candidates in a search for clips that contain people. Sound category and classifi- 
es cation scheme descriptors provide a means for organizing category concepts into hierarchical structures that enable 
this type of complex relational search strategy. 
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Sound Category 

[0035] As shown in Figure 6 for a simple taxonomy 600, a description sclieme (DS) is used for naming sound cate- 
gories. As an example, the sound of a dog barking can be given the qualitative category label "Dogs" 610 with "Bark" 
611 as a sub-category. In addition, "Woof" 612 or "Howl" 613 can be desirable sub-categories of "Dogs." The first two 
sub-categories are closely related, but the third is an entirely different sound event. Therefore, Figure 6 shows four 
categories are organized into a taxonomy with "Dogs" as the root node 61 0. Each category has at least one relation 
link 601 to another category in the taxonomy. By default, a contained category is considered a narrower category (NC) 
601 than the containing category. However, in this example, "Woof is defined as being a nearly synonymous with, but 
less preferable than, "Bark". To capture such structure, the following relations are defined as part of my description 
scheme. 

[0036] BC - Broader category means the related category is more general in meaning than the containing category. 

NC - Narrower category means the related category is more specific in meaning than the containing category. US - 
Use the related category that is substantially synonymous with the current category because it is preferred to the 
current category. UF - Use of the current category is preferred to the use of the nearly synonymous related category 
RC - The related category is not a synonym, quasi-synonym, broader or narrower category, but is associated with the 
containing category. 

[0037] The following XML-schema code shows how to instantiate the qualitative description scheme for the category 
taxonomy shown in Figure 6 using a description definition language (DDL): 



<SoundCategory tenn="r' scheme-"DOGS"> 
<Label>Dogs</Label> 

<TermRelation term="l.l" scheme="DOGS"> 
<Label>Bark</Label> 

<TennRelation term-" 1.2" scheme="DOGS" type="US"> 

<Label>Woof</Label> 
</TermRelation> 
</TermRelation> 

<TermRelation term-" 1.3" scheme="DOGS"> 

<Label>Howl</Label> 
</TermRelation> 
</SoxmdCategory> 



[0038] The category and scheme attributes together provide unique identifiers that can be used for referencing cat- 
egories and taxonomies from the quantitative description schemes, such as the probabilistic models described in great- 
er detail below. The label descriptor gives a meaningful semantic label for each category and the relation descriptor 
describes relationships amongst categorys in the taxonomy according to the invention. 

Classification Scheme 

[0039] As shown in Figure 7, categories can be combined by the relational links into a classification scheme 700 to 
make a richer taxonomy; for example, "Barks" 611 is a sub-category of "Dogs" 610 which is a sub-category of "Pets" 
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701 ; as is the category "Cats" 710. Cats 710 has the sound categories "IVIeow" 711 and "purr" 712. The following is an 
example of a simple classification scheme for "Pets" containing two categories: "Dogs" and "Cats". 
[0040] To implement this classification scheme by extending the previously defined scheme, a second scheme, 
named "CATS", is instantiated as follows: 



<SoundCategory tenTi="2" scheme="CATS"> 
<Label>Cats</Label> 

<TermRelation term="2.1 scheme="CATS"> 

<Label>Meow</Label> 
</TermRelation> 

<TennRelation term="2,2" scheme="CATS"> 

<Label>Purr</Label> 
</TerniRelatioii> 
</SouiidCategory> 



[0041] Now to combine these categories, a ClassificationScheme, called "PETS", is instantiated that references the 
previously defined schemes: 



<ClassificationScheme term- '0" scheme="PETS"> 
<Label>Pets</Label> 

<ClassificationSchemeRef scheme="DOGS"/> 
<ClassificationSchemeRef scheme="CATS"/> 
</ClassificationScheme> 



[0042] Now, the classifications scheme called "PETS" includes all of the category components of "DOGS" and "CATS" 
with the additional category "Pets" as the root. A qualitative taxonomy, as described above, is sufficient for text indexing 
applications. 

[0043] The following sections describe quantitative descriptors for classification and indexing that can be used to- 
gether with the qualitative descriptors to form a complete sound index and search engine. 

Quantitative Descriptors 

[0044] The sound recognition quantitative descriptors describe features of an audio signal to be used with statistical 
classifiers. The sound recognition quantitative descriptors can be used for general sound recognition including sound 
effects and musical instruments. In addition to the suggested descriptors, any other descriptor defined within the audio 
framework can be used for classification. 

Audio Spectrum Basis Features 

[0045] Among the most widely used features for sound classification are spectrum-based representations, such as 
power spectrum slices or frames. Typically, a each spectrum slice is an n-dimensional vector, with n being the number 
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of spectral channels, with up to 1 024 channels of data. A logarithmic frequency spectrum, as represented by an audio 
framework descriptor, helps to reduce the dimensionality to around 32 channels. Therefore, spectrum-derived features 
are generally incompatible with probability model classifiers due to their high dimensionality. Probability classifiers work 
best with fewer than 10 dimensions. 

[0046] Therefore, I prefer the low dimensionality basis functions produced by the single value decomposition (SVD) 
as described above and below. Then, an audio spectrum basis descriptor is a container for the basis functions that are 
used to project the spectrum to the lower-dimensional sub-space suitable for probability model classifiers. 
[0047] I determine a basis for each class of sound, and sub-classes. The basis captures statistically the most regular 
features of the sound feature space. Dimension reduction occurs by projection of spectrum vectors against a matrix 
of data-derived basis functions, as described above. The basis functions are stored in the columns of a matrix in which 
the number of rows corresponds to the length of the spectrum vector and the number of columns corresponds to the 
number of basis functions. Basis projection is the matrix product of the spectrum and the basis vectors. 

Spectrogram Reconstructed from Basis Functions 

[0048] Figure 8 shows a spectrogram 800 reconstructed from four basis functions according to the invention. The 
specific spectrogram is for "pop" music. The spectral basis vectors 801 on the left are combined with the basis projection 
vectors 802, using the vector outer product. Each resulting matrix of the outer product is summed to produce the final 
reconstruction. Basis functions are chosen to maximize the information in fewer dimensions than the original data. For 
example, basis functions can correspond to uncorrelated features extracted using principal component analysis (PGA) 
or a Karhunen-Loeve transform (KLT), or statistically independent components extracted by independent component 
analysis (ICA). The KLT or Hotelling transform is the preferred decorrelating transform when the second order statistics, 
i.e., covariances are known. This reconstruction is described in greater detail with reference to Figure 13. 
[0049] For classification purposes a basis is derived for an entire class. Thus, the classification space includes of 
the most statistically salient components of the class. The following DDL instantiation defines a basis projection matrix 
that reduces a series of 31 -channel logarithmic frequency spectra to five dimensions. 
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<AudioSpectrumBasis loEd.ge="62.5" hiEdge="8000" resolution="l/4 octave"> 


<Basis> 




<Matrixdim="31 5"> 


0 26 -0 05 


0 01 -0 70 0 44 


0.34 0.09 


0.21 -0.42 -0.05 


0.33 0.15 


0.24 -0.05 -0.39 


0.33 0.15 


0.24 -0.05 -0.39 


0.27 0.13 


0.16 0.24-0.04 


0.27 0.13 


0.16 0.24-0.04 


0.23 0.13 


0.09 0.27 0.24 


0.20 0.13 


0.04 0.22 0.40 


0.17 0.11 


0.01 0.14 0.37 


</Matrix> 


<;/Basis> 




</AudioSpectrumBasis> 



[0050] The loEdge, hiEdge and resolution attributes give lower and upper frequency bounds of the basis functions, 
and the spacing of the spectral channels in octave-band notation. In the classification framework according to the 
invention, the basis functions for an entire class of sound are stored along with a probability model for the class. 

Sound Recognition Features 

[0051] Features used for sound recognition can be collected into a single description scheme that can be used for 
a variety of different applications. The default audio spectrum projection descriptors perform well in classification of 
many sound types, for example, sounds taken from sound effect libraries and musical instrument sample disks. 
[0052] The base features are derived from an audio spectrum envelope extraction process as described above. The 
audio spectrum projection descriptor is a container for dimension-reduced features that are obtained by projection of 
a spectrum envelope against a set of basis functions, also described above. For example, the audio spectrum envelope 
is extracted by a sliding window FFT analysis, with a resampling to logarithmic spaced frequncy bands. In the preferred 
embodiment, the analysis frame period is 10ms. However, a sliding extraction window of 30ms duration is used with 
a Hamming window. The 30ms interval is chosen to provide enough spectral resolution to roughly resolve the 62.5 Hz- 
wide first channel of an octave-band spectrum. The size of the FFT analysis window is the next-larger power-of-two 
number of samples. This means for 30ms at 32kHz there are 960 samples but the FFT would be performed on 1024 
samples. For 30ms at 44.1 KHz, there are 1323 samples but the FFT would be performed on 2048 samples with out- 
of-window samples set to 0. 

[0053] Figures 9a and 9b show three spectral basis components 901-903 for a time index 910, and the resulting 
basis projections 911-913 with a frequency index 920 for a "laughter" spectrogram 1000 in Figures lOa-b. The format 
here is similar to those shown in Figures 4 and 5. Figure 10a shows a log scale spectrogram of laughter, and Figure 
10b a spectrogram reconstruction. Both figures plot the time and frequency indices on the x- and y-axes respectively. 
[0054] In addition to the base descriptors, a large sequence of alternative quantitative descriptors can be used to 
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define classifiers that use special properties of a sound class, such as the harmonic envelope and fundamental fre- 
quency features that are often used for musical instrument classification. 

[0055] One convenience of dimension reduction as done by my invention, is that any descriptor based on a scalable 
series can be appended to spectral descriptors with the same sampling rate. In addition, a suitable basis can be com- 
5 puted for the entire set of extended features in the same manner as a basis based on the spectrum. 

Spectrogram Summarization with a Basis Function 

[0056] Another application for the sound recognition features description scheme according to the invention is effi- 
10 cient spectrogram representation. For spectrogram visualization and summarization purposes, the audio spectrum 
basis projection and the audio spectrum basis features can be used as a very efficient storage mechanism. 
[0057] In order to reconstruct a spectrogram, we use Equation 2, described in detail below. Equation 2 constructs a 
two-dimensional spectrogram from the cross product of each basis function and its corresponding spectrogram basis 
projection, also as shown in Figure 8 as described above. 

15 

Probability IVIodel Description Schemes 
Finite State Model 

20 [0058] Sound phenomena are dynamic because spectral features vary over time. It is this very temporal variation 
that gives acoustic signals their characteristic "fingerprints" for recognition. Hence, my model partitions the acoustic 
signal generated by a particular source or sound class into a finite number of states. The partitioning is based on the 

spectral features. Individual sounds are described by their trajectories through this state space. This model is described 
in greater detail below with respect to Figures 11 a-b. Each state can be represented by a continuous probability 

25 distribution such as a Gaussian distribution. 

[0059] The dynamic behavior of a sound class through the state space is represented by a /c x /c transition matrix 
that describes the probability of transition to a next state given a current state. A transition matrix Tmodels the probability 
of transitioning from state /at time t-1 to state; at time An initial state distribution, which \sakx^ vector of probabilities, 
is also typically used in a finite-state model. The kth element in this vector is the probability of being in state k in the 

30 first observation frame. 

Gaussian Distribution Type 

[0060] A multi-dimensional Gaussian distribution is used for modeling states during sound classification. Gaussian 
35 distributions are parameterized by a 1 x n vector of means m, and an n x n covariance matrix, K, where n is the 
number of features in each observation vector. The expression for computation of probabilities for a particular vector 
X, given the Gaussian parameters is: 



40 



exp 



~(x-ni)^K-* (x-m) 



45 [0061] A continuous hidden Markov model is a finite state model with a continuous probability distribution model for 
the state observation probabilities. The following DDL instantiation is an example of the use of probability model de- 
scription schemes for representing a continuous hidden Markov model with Gaussian states. In this example, floating- 
point numbers have been rounded to two decimal places for display purposes only. 



55 
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<ProbabilityModel xsi:type="ContinuousMarkovModelType" numberStates="7"> 
<Imtial dim="7"> 
0.04 0.34 0.12 0.04 0.34 
<Transitions dim="7 7"> 
0.91 0.02 0.00 0.00 0.05 
0.01 0.99 0.00 0.00 0.00 
0.01 0.00 0.92 0.01 0.01 
0.00 0.00 0.00 0.99 0.01 
0.02 0.00 0.00 0.00 0.97 
0.00 0.00 0.01 0.00 0.00 
0.02 0.00 0.00 0.00 0.00 
</Transitions> 

<State><Label> 1</Label></State> 
<! — State 1 Observation Distribution ~ > 

<ObservationDistribution xsi;type=" GaussianDistributionType"> 
<Mean dim=="6"> 
5.11 -9.28 -0.69 -0.79 0.38 0.47 

</Mean> 

<Covariance dim— '6 6"> 
1.40-0.12-1.53-0.72 0.09-1.26 
-0.12 0.19 0.02-0.21 0.23 0.17 
-1.53 0.02 2.44 1.41 -0.30 1.69 
-0.72 -0.21 1.41 2.27-0.15 1.05 
0.09 0.23 -0.30-0.15 0.80 0.29 
-1.26 0.17 1.69 1.05 0.29 2.24 
</Covariance> 

<State><Label>2</Label></State> 
<! — Remaining states use same structures— > 
<\PobabilityModel> 



0.12 0.00 </Initial> 

0.01 0.01 

0.00 0.00 

0.06 0.00 

0.00 0.00 

0.00 0.00 

0.98 0.01 

0.02 0.96 



[0062] In this example, "ProbabilityModel" is instantiated as a Gaussian distribution type, which is derived from the 
base probability model class. 

Sound Recognition iUodei Description Schemes 

[0063] So far, I have isolated tools without any application structure. The following data types combine the above 

described descriptors and description schemes into a unified framework for sound classification and indexing. Sound 
segments can be indexed with a category label based on the output of a classifier. Additionally, the probability model 
parameters can be used for indexing sound in a database. Indexing by model parameters, such as states, is necessary 
for query-by-example applications when the query category is unknown, or when a narrower match criterion than the 
scope of a category is required. 

Sound Recognition iUiodei 

[0064] A sound recognition model description scheme specifies a probability model of a sound class, such as a 
hidden Markov model or Gaussian mixture model. The following example Is an instantiation of a hidden l\^arkov model 
of the "Barks" sound category 611 of Figure 6. A probability model and associated basis functions for the sound class 
is defined in the same manner as for the previous examples. 
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<SoundRecognitionModel id="sfxl.r' SoundCategoryRef^"Bark"> 
<ExtractionInfonnation tenn="Paraineters" scheme="ExtractionParameters"> 
<Label>NumStates=7, NumBasisComponents=5</Label> 
</ExtractionInfonTiation> 

<ProbabilityModel xsi:type="ContinuousMarkovModelType'* numberStates="7"> 

...<!— see previous example — > 
</ProbabilityModel> 

<SpectrumBasis loEdge="62.5" hiEdge="8000" resolution^" 1/4 octave"> 

...<!— see previous example — > 
</SpectruiiiBasis> 

</S ouxidRecQ gnitionMo del> 



Sound Model State Path 

20 

[0065] This descriptor refers to a finite-state probability model and describes the dynamic state path of a sound 
through the model. The sounds can be indexed in two ways, either by segmenting the sounds into model states, or by 
sampling of the state path at regular intervals. In the first case, each audio segment contains a reference to a state, 
and the duration of the segment indicates the duration of activation for the state. In the second case, the sound is 
25 described by a sampled series of indices that reference the model states. Sound categories with relatively long state- 
durations are efficiently described using the one-segment, one-state approach. Sounds with relatively short state du- 
rations are more efficiently described using the sampled series of state indices. 

[0066] Figure 11 a shows a log spectrogram (frequency v. time) 1100 of the dog-bark sound 611 of Figure 6. Figure 
1 1 b shows a sound model state path sequence of states through a continuous hidden Markov model for the bark model 
30 of Figure 11a, over the same time interval. In Figure lib, the x-axis is the time index, and the y-axis the state index. 

Sound Recognition Classifier 

[0067] Figure 12 shows a sound recognition classifier that uses a single database 1200 for all the necessary com- 
35 ponents of the classifier. The sound recognition classifier describes relationships between a number of probability 
models thus defining an ontology of classifiers. For example, a hierarchical recognizer can classify broad sound classes, 
such as animals, at the root nodes and finer classes, such as dogs:bark and cats:meow, at leaf nodes as described 
for Figures 6 and 7. This scheme defines mapping between an ontology of classifiers and a taxonomy of sound cate- 
gories using the graph's descriptor scheme structure to enable hierarchical sound models to be used for extracting 
40 category descriptions for a given taxonomy. 

[0068] Figure 13 shows a system 1300 for building a database of models. The system shown in Figure 13 is an 
extension of the system shown in Figure 1 . Here, the input acoustic signal is windowed before filtering to extract the 
spectrum envelope. The system can take audio input 1301 in the form of, e.g., WAV format audio files. The system 
extracts audio features from the files, and trains a hidden Markov model with these features. The system also uses a 
45 directory of sound examples for each sound class. The hierarchical directory structure defines an ontology that corre- 
sponds to a desired taxonomy. One hidden Markov model is trained for each of the directories in the ontology. 

Audio Feature Extraction 

50 [0069] The system 1300 of Figure 13 shows a method for extracting audio spectrum basis functions and features 
from an acoustic signal as described above. An input acoustic signal 1301 can either be generated by a single source, 
e.g., a human, or an animal, or a musical instrument, or a many sources, e.g., a human and an animal and multiple 
instruments, or even synthetic sounds. In the later case, the acoustic signal is a mixture. The input acoustic signal is 
first windowed 1310 into 1 0 ms frames. Note, in Figure 1 the input signal is band-pass filtered before windowing. Here, 

55 the acoustic signal is first windowed and than filtered 1320 to extract a short-time logarithmic-in-frequency spectrum. 
The filtering performs a time-frequency power spectrum analysis, such as a squared-magnitude short-time Fourier 
transform. The result is a matrix with A/7 frames and A/frequency bins. The spectral vectors x, are the rows of this matrix. 
[0070] Step 1330 performs log-scale normalization. Each spectral vector x is converted from the power spectrum to 
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a decibel scale 1331 z = 10 log (x). Step 1332 determined the L2-norm of the vector elements 



. The new unit-norm spectral vector is then determined the spectrum envelope^ by z/r, which divides each slice z by 
its power r, and the resulting normalized spectrum envelope ^ 1340 is passed to the basis extraction process 1360. 
[0071] The spectrum envelope ^ places each vector row-wise in the form of an observation matrix. The size of the 
resulting matrix is MxA/ where /Wis the number 

of time frames and A/ is the number of frequency bins. The matrix will have the following structure: 




Basis Extraction 

[0072] Basis functions are extracted using the singular value decomposition SVD 130 of Figure 1 . The SVD is per- 
formed using the command [U, S, V] = SVD(X, 0). I prefer to use an "economy" SVD. An economy SVD omits unnec- 
essary rows and columns during the factorization of the SVD. I do not need the row-basis functions, thus the extraction 
efficiency of the SVD is increased. The SVD factors the matrix as follows. ^ = USV"^, where ^ is factored into a matrix 
product of three matrices, the row basis U, the diagonal singular value matrix S, and the transposed column basis 
functions V. The basis is reduced by retaining only the first /C basis functions, i.e., the first columns of V: 

where K \s typically in the range of 3-10 basis functions for sound feature-based applications. To determine the pro- 
portion of information retained for K basis functions use the singular values contained in matrix S: 



where l(K) is the proportion of information retained for /< basis functions, and N is the total number of basis functions 
which is also equal to the number of spectral bins. The SVD basis functions are stored in the columns of the matrix. 
[0073] For maximum compatibility between applications, the basis functions have columns with unit L2-norm, and 
the functions maximize the information in /c dimensions with respect to other possible basis functions. Basis functions 
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can be orthogonal, as given by PGA extraction, or non-ortinogonal as given by ICA extraction, see below. Basis pro- 
jection and reconstruction are described by the following analysis-synthesis equations, 

5 Y = XV (1) 

and 

10 X = YV", (2) 

where X is the spectrum envelope, Y are the spectral features, and V are the temporal features. The spectral features 

are extracted from the m x /(observation matrix of features, X is the m x n spectrum data matrix with spectral vectors 
organized row wise, and V is a r? x /( matrix of basis functions arranged in the columns. 

[0074] The first equation corresponds to feature extraction and the second equation corresponds to spectrum re- 
construction, see Figure 8, where V+ denotes the pseudo inverse of V for the non-orthogonal case. 

Independent Component Analysis 

[0075] After the reduced SVD basis V has been extracted, an optional step can perform a basis rotation to directions 

of maximal statistical independence. This isolates independent components of a spectrogram, and is useful for any 
application that requires maximum separation of features. To find a statistically independent basis using the basis 
functions obtained above, any one of the well-known, widely published independent component analysis (ICA) proc- 
esses can be used, for example, JADE, or FastlCA, see Cardoso, J.F. and Laheld, B.H. "Equivariant adaptive source 
separation," IEEE Trans. On Signal Processing, 4:112-114, 1996, or Hyvarinen, A. "Fast and robust fixed-point algo- 
rithms for independent component analysis," IEEE Trans. On Neural Networks, 1 0(3):626-634, 1999. 
[0076] The following use of ICA factors a set of vectors into statistically independent vectors [V^,A]=ica(V^;f), where 
the new basis is obtained as the product of the SVD input vectors and the pseudo-inverse of tlie estimated mixing 
matrix A given by the ICA process. The ICA basis is the same size as the SVD basis and is stored in the columns of 
the basis matrix. The retained information ratio, l(K), is equivalent to the SVD when using the given extraction method. 
The basis functions V^ 1 361 can be stored in the data base 1 200. 

[0077] In the case where the input acoustic signal is a mixture generated from multiple sources, the set of features 
produced by the SVD can be clustered into groups using any known clustering technique having a dimensionality equal 
to the dimensionality of the features. This puts like features into the same group. Thus, each group includes features 

for the acoustic signal generate by a single source. 

[0078] The number of groups to be used in the clustering can be set manually or automatically, depending on a 
desired level of discrimination desired. 

Use of Spectrum Subspace Basis Functions 

40 

[0079] To obtain projection or temporal features Y, the spectrum envelope matrix X is multiplied by the basis vectors 
of the spectral features V. This step is the same for both for SVD and ICA basis functions, i.e., =xv, where Y is a 
matrix consisting of the reduced dimension features after projection of the spectrum against the basis V. 
[0080] For independent spectrogram reconstruction and viewing, I extract the non-normalized spectrum projection 
by skipping the normalization step 1330 extraction, thus, Y;^ = XV;^. Now, to reconstruct an independent spectrogram, 
Xk as shown in Figure 8, component use the individual vector pairs, corresponding to the Kth projection vector and 
the inverted Kth basis vector and apply the reconstruction equation Xi^ = y^ where the operator indicates the 
transpose for SVD basis functions, which are orthonormal, or the pseudo-inverse for ICA basis functions, which is non- 
orthogonal. 

50 

Spectrogram Summarization by Independent Components 

[0081] One of the uses for these descriptors is to efficiently represent a spectrogram with much less data than a full 
spectrogram. Using an independent component basis, individual spectrogram reconstructions, e.g., as seen in Figure 
8, generally correspond to source objects in the spectrogram. 



14 



EP 1 260 968 A1 



Model Acquisition and Training 

[0082] Much of the effort in designing a sound classifier is spent collecting and preparing training data. The range 
of sounds should reflect the scope of the sound category. For example, dog barks can include individual barks, multiple 
5 barks in succession, or many dogs barking at once. The model extraction process adapts to the scope of the data, 
thus a narrower range of examples produces a more specialized classifier. 

[0083] Figure 14 show a process 1400 for extracting features 1410 and basis function 1420, as described above, 
from acoustic signals generated by known sources 1401. These are then used to train 1440 hidden Markov models. 
The trained models are stored in the database 1200 along with their corresponding features. During training, an un- 
10 supervised clustering process is used to partition an n-dimensional feature space into /estates. The feature space is 
populated by reduced-dimension observation vectors. The process determines an optimal number of states for the 
given data by pruning a transition matrix given an initial guess for k. Typically, between five and ten states are sufficient 
for good classifier performance. 

[0084] The hidden Markov models can be trained with a variant of the well-known Baum-Welch process, also known 
15 as Fonward- Backward process. These processes are extended by use of an entropic prior and a deterministic annealing 
implementation of an expectation maximization (EM) process. 

[0085] Details for a suitable HMM training process 1430 are described by Brand in "Pattern discovery via entropy 
minimization" In Proceedings, Uncertainty'99. 

[0086] Society of Artificial intelligence and Statistics #7, Morgan Kaufmann, 1999, and Brand, "Structure discovery 
20 in conditional probability models via an entropic prior and parameter extinction, " Neural Computation, 1 999. 

[0087] After each HMM for each known source is trained, the model is saved in permanent storage 1200, along with 
its basis functions, i.e., the set of sound features. When a number of sound models have been trained, corresponding 
to an entire taxonomy of sound categories, the HMMs are collected together into a larger sound recognition classifier 
data structure thereby generating an ontology of models as shown in Figure 12. The ontology is used to index new 
25 sounds with qualitative and quantitative descriptors. 

Sound Description 

[0088] Figure 15 shows an automatic extraction system 1500 for indexing sound in a database using pre-trained 
30 classifiers saved as DDL files. An unknown sound is read from a media source format, such as a WAV file 1501 . The 
unknown sound is spectrum projected 1 520 as described above. The projection, that is, the set of features is then used 
to select 1530 one of the HMMs from the database 1200. A Viterbi decoder 1540 can be used to give both a best-fit 
model and a state path through the model for the unknown sound. That is, there is one model state for each windowed 
frame of the sound, see Figure 11b. Each sound is then indexed by its category, model reference and the model state 
35 path and the descriptors are written to a database in DDL format. The indexed database 1599 can then be searched 
to find matching sounds using any of the stored descriptors as described above, for example, all dog barkings. The 
substantially similar sounds can then be presented in a result list 1560. 

[0089] Figure 16 shows classification performance for ten sound classes 1601-1610, respectively: bird chirps, ap- 
plause, dog barks, explosions, footsteps, glass breaking, gun shots, gym shoes, laughter, and telephones. Perform- 
40 ance of the system was measured against a ground truth using the label of the source sound as specified by a pro- 
fessional sound-effect library. The results shown are for novel sounds not used during the training of the classifiers, 
and therefore demonstrate the generalization capabilities of the classifier. The average performance is about 95% 
correct. 

"^5 Example Search Applications 

[0090] The following sections give examples of how to use the description schemes to perform searches using both 
DDL-based queries and media source-format queries. 

50 Query by Example with DDL 

[0091] As shown in Figure 17 in simplified form, a sound query is presented to the system 1700 using the sound 
model state path description 1710 in DDL format. The system reads the query and populates internal data structures 
with the description information. This description is matched 1550 to descriptions taken from the sound database 1599 
55 stored on disk. The sorted result list 1560 of closest matches is returned. 

[0092] The matching step 1 550 can use the sum of square errors (SSE) between state-path histograms. This match- 
ing procedure requires little computation and can be computed directly from the stored state-path descriptors. 
[0093] State-path histograms are the total length of time a sound spends in each state divided by the total length of 
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the sound, thus giving a discrete probability density function with the state index as the random variable. The SSE 
between the query sound histogram and that of each sound in the database is used as a distance metric. A distance 
of zero implies an identical match and increased non-zero distances are more dissimilar matches. This distance metric 
is used to rank the sounds in the database in order of similarity, then the desired number of matches is returned, with 

the closest match listed first. 

[0094] Figure 1 8a shows a state path, and Figure 1 8b a state path histogram for a laughter sound query. Figure 1 9a 
shows state paths and Figure 1 9b histograms for the five best matches to the query. All matches are from the same 
class as the query which indicates the success the correct performance of the system. 

[0095] To leverage the structure of the ontology, sounds within equivalent or narrower categories, as defined by a 
taxonomy, are returned as matches. Thus, the 'Dogs' category will return sounds belonging to all categories related to 
'Dogs' in a taxonomy. 

Query-by-Example with Audio 

[0096] The system can also perform a query with an audio signal as input. Here, the input to the query-by-example 
application is an audio query instead of a DDL description-based query. In this case, the audio feature extraction 
process is first performed, namely spectrogram and envelope extraction is followed by projection against a stored set 
of basis functions for each model in the classifier. 

[0097] The resulting dimension-reduced features are passed to the Viterbi decoder for the given classifier, and the 
HMM with the maximum-likelihood score for the given features is selected. The Viterbi decoder essentially functions 
as a model-matching algorithm for the classification scheme. The model reference and state path are recorded and 
the results are matched against a pre-computed database as in the first example. 

[0098] It is to be understood that various other adaptations and modifications may be made within the spirit and 
scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications 
as come within the true spirit and scope of the invention. 



Claims 

1. A method for extracting features from an acoustic signal generated from a single source, comprising: 

windowing and filtering the acoustic signal to produce a spectral envelope; and 

reducing the dimensionality of the spectral envelope to produce a set of features, the set including spectral 
features and corresponding temporal features characterizing the single source. 

2. The method of claim 1 further comprising: 

multiplying the spectral features and temporal features using a outer product to reconstruct a spectrogram of 
the accoustic signal. 

3. The method of claim 1 further comprising: 

applying independent component analysis to the set of feature to separate the features in the set. 

4. The method of claim 1 further comprising: 

log-scaling and L2-normalizing the spectral envelope to a decibel scale and unit L2-norm before reducing the 

dimensionality of the spectral envelope. 

5. A method for extracting features from an acoustic signal generated from a plurality of sources, comprising: 

windowing and filtering the acoustic signal to produce a spectral envelope; 
reducing the dimensionality of the spectral envelope to produce a set of features; 

clustering the features in the set to produce a group of features for each of the plurality of sources, the features 
in each group including spectral features and corresponding temporal features characterizing each source. 

6. The method of claim 5 wherein each group of features is a quantitative descriptor of each source, and futher 
comprising: 
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associating a qualitative descriptor witin eacli quantitative descriptor to generate a category for each source. 



7. The method of claim 6 further comprising: 



5 organizing the categories in a database as a taxonomy of classified sources; 

relating each category with at least one other category in the database by a relational link. 

8. The method of claim 7 wherein the categories are stored in the database using a description definition language. 

10 9. The method of claim 8 wherein a particular category in a DDL instantiation defines a basis projection matrix that 
reduces a series of logarithmic frequencies spectra of a particular source to fewer dimensions. 

10. The method of claim 6 wherein the categories include environmental sounds, bacl<ground noises, sound effects, 
sound textures, animal sounds, speech, non-speech utterances, and music. 

15 

1 1 . The method of claim 7 further comprising: 

combining substantially similar categories in the database as a hierarchy of classes. 

20 12. The method of claim 6 a particular quantitative descriptor further includes a harmonic envelope descriptor, and 
fundamental frequency descriptor. 



13. The method of claim 5 wherein the temporal features describe a trajectory of the spectral features over time, and 
further comprising: 

25 

partitions the acoustic signal generated by a particular source into a finite number of states based on the 

corresponding spectral features; 

representing each state by a continuous probability distribution; 

representing the temporal features by a transition matrix to model probabilities of transitions to a next state 
30 given a current state. 



14. The method of claim 13 wherein the continuous probability distribution is a Gaussian distribution parameterized 
by a 1 X n vector of means m, and an n x n co variance matrix K, where n is the number of spectral features in 
each spectral envelope, and the probabilities of a particular spectral envelope x is given by: 

35 



1 



-exp 



(27tY\K\^ 



in)^K-^(x-ni) 



40 



15. The method of claim 5 wherein each source is known, and further comprising: 



training, for each known source, a hidden Markov model with the set of features; 
45 storing each trained hidden Markov model with the associated set of spectral features in a database. 

16. The method of claim 5 wherein a set of acoustic signals belongs to a known category, and further comprising: 

extracting a spectral basis for the acoustic signals; 
50 training a hidden Markov model using the temporal features of the acoustic signals; 

storing each trained hidden Markov model with the associated spectral basis features. 



17. The method of claim 15 further comprising: 



55 generating an unknown acoustic from an unknown source; 

windowing and filtering the unknown acoustic signal to produce an unknown spectral envelope; 

reducing the dimensionality of the unknown spectral envelope to produce a set of unknown features, the set 

including unknown spectral features and corresponding unknown temporal features characterizing the un- 



17 



EP 1 260 968 A1 



known source; 

selecting one of the stored hidden Markov models that best-fits the unknown set of features to identify the 
unknown source. 

18. The method of claim 1 7 wherein a plurality of the stored hidden IVIarkov models are selected to identify a plurality 
of known source substantially similar to the unknown source. 
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